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ABSTRACT 

We combine previously published interferometric and single-dish data of relatively 
nearby massive dense cores that are actively forming stars to test whether their ‘frag¬ 
mentation level’ is controlled by turbulent or thermal support. We hnd no clear correla¬ 
tion between the fragmentation level and velocity dispersion, nor between the observed 
number of fragments and the number of fragments expected when the gravitationally 
unstable mass is calculated including various prescriptions for ‘turbulent support’. On 
the other hand, the best correlation is found for the case of pure thermal Jeans frag¬ 
mentation, for which we infer a core formation efficiency around 13%, consistent with 
previous works. We conclude that the dominant factor determining the fragmentation 
level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans 
fragmentation. 

Key words: stars: formation, clusters — ISM: lines and bands — ISM: structure — 
radio continuum: ISM — turbulence 


et al. 2014; Federrath 2015; Guszejnov & Hopkins 2015; 
Salim, Federrath, & Kewley 2015). However, some recent 
studies suggest that this may not be the case. For instance, 
molecular clouds seem to form by large-scale compressions 
in the diffuse, warm, HI medium. The compressed gas un¬ 
dergoes a transition to the cold, dense atomic phase (e. g., 
Hennebelle & Perault 1999, Heitsch et al. 2005, Vazquez- 
Semadeni et al. 2006), which is highly prone to Jeans in¬ 
stability (Hartmann et al. 2001), and thus must begin soon 
to collapse, in spite of the turbulence generated inside it by 
the original compression (Koyama & Inutsuka 2002, Audit & 
Hennebelle 2005, Vazquez-Semadeni et al. 2007, Heitsch & 
Hartmann 2008). Moreover, once molecular clouds achieve 
column densities ~ 10^^ cm“^, they are able to form molecu¬ 
lar gas (Bergin et al. 2004) so that the formation of molecules 
may be essentially a byproduct of the gravitational col¬ 
lapse of the clouds (Hartmann et al. 2001). In this alterna¬ 
tive scenario, the observed non-thermal motions of molec- 


1 INTRODUCTION 

For more than 60 years it has been thought that turbu¬ 
lence is an agent capable of providing support to molecular 
clouds against gravitational collapse (e. g., Chandrasekhar 
1951, Bonazzola et al. 1987, Leorat et al. 1990, McKee & 
Tan 2003), while simultaneously producing local density en¬ 
hancements that may become Jeans-unstable and collapse 
(e. g., Sasao 1973, Elmegreen 1993, Padoan 1995, Vazquez- 
Semadeni & Gazol 1995, Klessen et al. 2000, Vazquez- 
Semadeni et al. 2003), and this is currently the most ac¬ 
cepted scenario for the dynamical state of molecular clouds 
(e. g., Mac Low & Klessen 2004; Krumholz & McKee 2005; 
Hennebelle & Chabrier 2008, 2011; Hopkins 2012; Chabrier 
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ular clouds, rather than consisting of random, small-scale, 
isotropic motions that can act as a pressure, would actually 
be dominated by inward motions caused by the gravitational 
collapse, which would occur both at large and small scales 
in a hierarchical and chaotic fashion (Vazquez-Semadeni et 
al. 2009, Ballesteros-Paredes et al. 2011a). This implies that 
the bulk of the observed non-thermal motions cannot pro¬ 
vide support against the self-gravity of the clouds. 

In previous papers of this series, we have presented ev¬ 
idence that the dynamics of molecular clouds are indeed 
dominated by gravity by showing that this scenario unihes 
molecular clouds and massive clumps in a single scaling rela¬ 
tion (Heyer et al. 2009) that extends those by Larson (1981, 
Ballesteros-Paredes et al. 2011a), and by showing that nu¬ 
merical simulations of cloud evolution including self-gravity 
develop power-law high-density tails in their column density 
probability distribution functions as a consequence of the 
gravitational collapse (see also Klessen 2000, Kritsuk et al. 
2011, Ballesteros-Paredes et al. 2011b; Federrath & Klessen 
2013), in agreement with observations (e. g., Kainulainen et 
al. 2009; Schneider et al. 2013). In the present contribution, 
we present a further line of evidence, by examining the mech¬ 
anism responsible for fragmentation of dense cores. Indeed, 
a still unsolved and highly debated question is what are the 
main drivers of fragmentation in massive dense core^, which 
are believed to be the precursors of stellar clusters. The cru¬ 
cial parameter to estimate the fragmentation level of a dense 
core is the Jeans mass, which in its general form takes into 
account the different mechanisms of support against grav¬ 
ity, through the use of the ‘effective sound speed’, Cefi (e. g., 
Mac Low & Klessen 2004). Among the most debated forms 
of support are turbulent and thermal support. Thus, if the 
average turbulence level, average temperature and average 
density of a massive dense core are known, one can easily cal¬ 
culate the Jeans mass in both cases and estimate, given the 
mass of the core, the number of fragments expected in each 
case, so that one can assess which form of support against 
gravity is controlling the fragmentation process. 

Up to now, this simple question could not be answered 
because of a lack of statistically significant samples of mas¬ 
sive dense cores where the fragmentation level has been 
assessed in a uniform way and down to spatial resolu¬ 
tions comparable to separations between cluster members 
(~ 1000 AU). Recently, Palau et al. (2013, 2014) compiled 
a sample of 19 massive dense cores with on-going star for¬ 
mation and studied the fragmentation level within the cores 
down to ~ 1000 an, and ~ 0.5 Mq of mass sensitivity. This 
would be a first approach to study the number of protostars 
(compact fragments will most likely become protostars, see 
Palau et al. 2013, 2014) within a massive dense core (as¬ 
sumed to become a cluster). In these works, the fragmen¬ 
tation level was assessed by counting the number of mil¬ 
limetre sources within a field of view of 0.1 pc of diameteiQ, 

^ We will follow the nomenclature of Williams et al. (2000) and 
Bontemps et al. (2010) where a massive dense core refers to a 
dense gas structure of ~ 0.1 pc in size and > 20 Mq in mass, 
which does not necessarily collapse into one star but can fragment 
into compact condensations and form a small cluster of stars. 

^ The size of 0.1 pc was taken to be the smallest size where 
fragmentation could be studied (because of the limitations by the 
primary beam response) in the works of Palau et al. (2013, 2014). 


Amm. Furthermore, by fitting the spectral energy distribu¬ 
tions and submillimetre intensity profiles of the cores, Palau 
et al. (2014) modelled their density and temperature struc¬ 
ture, so that densities and temperatures at different spatial 
scales could be estimated and thermal Jeans masses could 
be calculated. In this work, we compile observational data 
based on dense gas tracers for the sample of cores of Palau 
et al. (2013, 2014) and analyze them in a uniform way, in 
order to assess the turbulence level, estimate the turbulent 
Jeans mass, and Hnally compare the fragmentation level ob¬ 
served to the fragmentation level expected for each form of 
support, turbulent or thermal. 

The plan of the present contribution is as follows: 321 
presents the compiled data we used. In 33] we present the 
fragmentation level of our cores, and how this fragmentation 
correlates (or not) with their physical parameters. Finally, 
in 331 we discuss the physical implications of our results and 
present our main conclusions. 


2 THE SAMPLE AND DATA COMPILATION 

The present work is based on the sample of massive dense 
cores presented in Palau et al. (2013, 2014), whose distances, 
luminosities, and masses range from 0.45 to 3 kpc, from 300 
to 10® Z/0, and from 80 to 1500 Mq, respectively (given 
in Table Al). The sample was selected from deeply embed¬ 
ded intermediate/high-mass star-forming regions published 
in the literature that have been studied in the millimetre 
range down to mass sensitivities of ~ 0.5 Mq, and spatial 
resolutions of ~ 1000 au. Palau et al. (2014) modelled the 
massive dense cores with temperature and density proHles 
decreasing with radius following power-laws, and determined 
a number of properties of the density and temperature struc¬ 
ture in a uniform way for all the sample. 

In order to estimate the turbulence level in each core 
in a uniform way and compare it to the ‘turbulent fragmen¬ 
tation’ level, we used Very Large Array (VLA) NH3(1,1) 
data in C/D conhguration, available for 14 (out of 19) mas¬ 
sive dense cores (Torrelles et al. 1989; Zhou et al. 1990; 
Mangum et al. 1992; Tieftrunk et al. 1998; Wiseman & Ho 
1998; Gomez et al. 2003; Sanchez-Monge et al. 2013; see 
also Palau et al. 2014). VLA beams are typically < 5 arcsec, 
and the minimum angular scales filtered out by the inter¬ 
ferometer are > 35 arcsec (estimated following appendix in 
Palau et al. 2010, and using a minimum baseline of 35 m). 
The latter corresponds to ~ 0.3 pc (for typical distances of 
the regions in the sample), slightly larger than the field of 
view of 0.1 pc that we are studying. Thus, with these VLA 
NHa data we are recovering most of the emission at the spa¬ 
tial scales we are studying, with an angular resolution good 
enough (~ 5 arcsec) to resolve typical sizes of massive dense 
cores (e.g., Sanchez-Monge et al. 2013). The NH3(1,1) hy- 
perfine structure was fitted and hyperfine ‘observed’ FWHM 
line widths, Auobs, were inferred for each core. We note that 
we took special care to make the sample as uniform as pos¬ 
sible and we measured the NH3(1,1) line width in all cases 
by using the ‘nh3(l,l) method’ in the CLASS package of the 
GILDAS software, on the spectrum resulting from averaging 
the NH3(1,1) emission over the central region of ~ 0.1 pc of 
diameter of the dense core where we study the fragmentation 
level. When literature did not provide an average spectrum 
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Table 1. Modelled properties of the massive dense cores and compiled observational velocity dispersions 


Source 

Nmm 

-^O.lpc^ 

(Mq) 

(10^ cm“^) 

To.ipc'^ 

(K) 

NH3 b 
^ lD,obs 
(km s“^) 

NH3 b 
^ ID,nth 
(km s“^) 

A4nh3*’ 

N2H+ b 
^lD,obs 
(km s“^) 

N2H+ b 
^ ID,nth 
(km s“^) 

AIN2H+*’ 

1-IG1396N 

4 

11 

3.6 

25 

_ 

_ 

_ 

0.79 

0.78 

4.5 

2-122198== 

1.5 

11 

3.6 

26 

0.47 

0.45 

2.6 

0.59 

0.59 

3.4 

3-N2071-1 

4 

17 

5.7 

24 

0.44 

0.43 

2.5 

- 

- 

- 

4-N7129-2 

1 

11 

3.6 

35 

- 

- 

- 

0.59 

0.58 

2.8 

5-CB3-mm 

2 

15 

5.2 

40 

- 

- 

- 

0.73 

0.72 

3.3 

6-I22172N 

3 

9 

3.2 

48 

0.59 

0.58 

2.4 

0.87 

0.86 

3.6 

7-OMC-lS 

9 

38 

13 

49 

1.11 

1.10 

4.5 

0.90 

0.89 

3.7 

8-A5142 

7 

39 

13 

47 

1.61 

1.61 

6.8 

1.09 

1.08 

4.6 

9-I05358NE 

4 

27 

9.1 

35 

0.72 

0.71 

3.5 

1.07 

1.07 

5.3 

10-120126 

1 

14 

4.8 

68 

2.00 

1.99 

7.0 

0.85 

0.84 

2.9 

11-122134 

3.5 

10 

3.2 

50 

0.71 

0.70 

2.8 

0.62 

0.61 

2.5 

12-HH80-81 

3 

12 

4.2 

66 

0.74 

0.72 

2.6 

- 

- 

- 

13-W3IRS5 

3.5 

12 

4.0 

138 

0.87 

0.84 

2.1 

1.18 

1.17 

2.9 

14-A2591 

1.5 

16 

5.2 

147 

0.68 

0.62 

1.5 

- 

- 

- 

15-Gyg-N53 

6 

30 

10 

27 

0.17=* 

0.13=^ 

0.7=^ 

0.81 

0.80 

4.4 

16-Gyg-N12 

2.5 

15 

5.0 

29 

- 

- 

- 

1.23 

1.23 

6.6 

17-Gyg-N63 

1 

14 

4.6 

31 

- 

- 

- 

0.82 

0.82 

4.2 

18-Gyg-N48 

5 

35 

12 

36 

1.25 

1.25 

6.1 

1.21 

1.21 

5.9 

19-DR21-OH 

11 

69 

23 

49 

1.51 

1.51 

6.3 

- 

- 

- 


^ A^O.lpc is the mass inside a region of 0.1 pc of diameter computed according: Mo.ipc = M{R = 0.05pc) = 47rpor‘Q ^ , where p, ro, 
and po are index of the density power law, the reference radius adopted to be 1000 AU, and the density at the reference radius (given 
in Table 1 of Palau et al. 2014); no.ipc To.ipc correspond to average density and temperature inside a region of 0.1 pc of diameter. 

f T(r)p(r)r^dr 

To.i pc is estimated as Tr = ° -, where T{r) and p(r) were calculated as power laws of the form T{r) = To{r/ro)~^ and 

J p(r)r^dr 

p(r) = po('^/^o)~^) with To and po being the values at the reference radius ro of 1000 AU. To, po, p, and q are given in Table 1 of Palau 
et al. (2014). The final expression is Tr = (l^) ^' 

'^ll?oba '^iD^ba Calculated from the measured FWHM line width, Auoba) 3-s = AuqUb/CS cm = 

\/°^lD,oba ~ °^th’ ™ith (7th = 's/T/ mn) (^B being the Boltzmann constant, fi the molecular weight (17 for NH3, 29 for N2H'^), 
rrtH the mass of the hydrogen atom and T the temperature of the region, taken from column (5) of this table). The Mach number M is 
calculated as (rsD.nth/ca, with Ca being the sound speed calculated as Cg = T/{fi mu), using fi = 2.3, and crsD.nth = y/S 

The parameters of the density and temperature structure for this source are different from Palau et al. (2014) because here we have 
used the original JCMT data of Jenness et al. (1995) and we have not assumed any error beam in the modelling (see Appendix B). 

Marginal detection of the NH3(1,1) line, not taken into account in the analysis of this work. 

Refs: IC1396N: Alonso-Albi et al. (2010); 122198: Sanchez-Monge et al. (2013), Fontani et al. (2011); NGC2071-1: Zhou et al. (1990); 
NGC7129-2: Fuente et al. (2005); CB3-mm: Alonso-Albi et al. (2010); I22172N: Sanchez-Monge et al. (2013); Fontani et al. (2006); OMG- 
IS: Wiseman &; Ho (1998); Tatematsu et al. (2008); A5142: Sanchez-Monge et al. (2013); Fontani et al. (2011); I05358NE: Sanchez-Monge 
et al. (2013); Fontani et al. (2011); 120126: Sanchez-Monge et al. (2013); Fontani et al. (2006); 122134: Sanchez-Monge et al. (2013); 
Fontani et al. (2015); HH80-81: Gomez et al. (2003); W3IRS5: Tieftrunk et al. (1998); Gerner et al. (2014); A2591: Torrefies et al. (1989); 
Cyg-N53: VLA archive; Bontemps et al. (2010); Gyg-N12: Bontemps et al. (2010); Cyg-N63: Bontemps et al. (2010); Cyg-N48: Mangum 
et al. (1992), Bontemps et al. (2010); DR21-OH: Mangum et al. (1992). 


over the massive dense core that we are studying, we down¬ 
loaded and reduced the VLA data to extract the spectrum. 
This was done for HH 80-81, W3IRS5, A2591, Cyg-N53, 
Cyg-N48, and DR21-OH (VLA projects AG0552, AT0180, 
AT0084, AW0240, AF386, respectively: standard calibration 
as described in Sanchez-Monge et al. 2013 was applied). In 
two sources (122134 and DR21-OH) we fitted two velocity 
components. From these line widths we calculated the ob¬ 
served velocity dispersions as criD,oba = Auoba/(8ln2)^/^. 
The values of the compiled observed velocity dispersions 
from NH3(1,1) VLA data are listed in Table[T] and the spec¬ 
tra with the corresponding fits to the hyperfine structure are 
shown in Appendix B. 


Because in some cases the VLA NH3(1,1) emission 
might be affected by the passage of an outflow (e. g., the 
NH3(1,1) line width of IRAS 20126-1-4104 is larger along the 
direction of the outflow, see Fig. B1 of Sanchez-Monge et 
al. 2013), we additionally compiled data from a different 
dense gas tracer, N 2 H"''(l- 0 ), observed using a single-dish 
telescope (IRAMSOm in all cases except for OMC-IS, for 
which Nobeyama45m was used). This is a reasonable ap¬ 
proach to avoid contamination by outflow for several rea¬ 
sons. First, N 2 H^ is known to be destroyed by CO (e. g., 
Joergensen 2004; Busquet et al. 2011). Second, the outflow is 
typically compact and thus its emission should be diluted in 
the single-dish beam. Thus, the N 2 H''" line widths should be 
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less affected by the passage of the outflow and thus more reli¬ 
able to measure the ‘initiar non-thermal motions unaffected 
by stellar feedback. The N 2 H’^(l- 0 ) data were available for 
15 out of 19 regions and were compiled from the literature 
(Fuente et al. 2005; Fontani et al. 2006, 2011, 2015; Tatem- 
atsu et al. 2008; Alonso-Albi et al. 2010; Bontemps et al. 
2010; Gerner et al. 2014), and its hyperfine structure was 
fitted using the CLASS package of the GILDAS software. 
The IRAMSOm Telescope provides a beam of ~ 26 arcsec 
at the frequency of N 2 H+(l- 0 ), comparable to the spatial 
scale at which the massive dense cores are being studied 
( 0.1 pc, at the typical distances of the cores of our sample), 
and about a factor of 5 larger than the VLA beam. By using 
the same method outlined above for NH 3 ( 1 , 1 ) we inferred 
the observed velocity dispersions for N 2 H"*’(l- 0 ), listed in 
Table [T] and the spectra and fits are shown in Appendix 
B. The method used to fit the hyperfine structure for both 
N2H''’ and NH3 takes the opacity effects into account. The 
observed N 2 H^(l- 0 ) velocity dispersions range from 0.6 to 
1.2 km s“^, a narrower range than that of the velocity dis¬ 
persions inferred from VLA NH3 data (ranging from 0.5 to 
2.0 km s“^), as expected (because of the outflow contami¬ 
nation of the NH3 VLA data). 

The main difference between the line widths reported 
in this work and the line widths reported in column ( 10 ) of 
Table 2 of Palau et al. (20l43) is that we here re-reduced the 
interferometric data and re-did the fits of the spectra in all 
cases using the same method, instead of just taking the val¬ 
ues reported in the literature, which use different methods. 
Thus the present analysis is uniform in the sense that the 
method to infer the line widths is the same for all sources. 


3 RESULTS AND ANALYSIS 

3.1 Fragmentation level vs. density and velocity 
dispersions 

The compiled observed velocity dispersions of NH 3 ( 1 , 1 ) and 
N 2 H"*’( 1 - 0 ) together with the modeling of the temperature 
structure of the massive dense cores (Palau et al. 2014) 
allowed us to separate the thermal from the non-thermal 
contribution of the observed velocity dispersion. We esti¬ 
mated the thermal comp onent of the velocity dispersion, 
ath, from ks T/{jj, mn), with ks the Boltzmann constant, 
jj. the molecular weight (17 for NH 3 , 29 for N 2 LL^), mn the 
mass of the hydrogen atom, and T the temperature of the re¬ 
gion, which was adopted from the average density-weighted 
temperature inside a region of 0.1 pc of diameter (the same 
region where we assessed the fragmentation level). This aver¬ 
age temperature is estimated from the density and tempera¬ 
ture power-laws modelled by Palau et al. (2014, see notes of 
Table [1] for further details) for each core. The non-thermal 
component was estimated by using uiD.nth = V'^obs “ ^th- 


® Linewidths reported in column (9) of Table 2 of Palau et al. 
(2014), or in column (10) of Table 4 of Palau et al. (2013) corre¬ 
spond to quiescent cores in the surroundings of the massive dense 
cores where fragmentation is being studied, and are not compara¬ 
ble to the line widths reported here, corresponding to the massive 
dense cores where active star formation is taking place and where 
fragmentation is being studied. 




‘^nth-NH3(l<m/s) 



Figure 1. Observed ‘fragmentation level’ (A^mm) vs different 
quantities (Table Q . a) A^mm vs the density of the core within a 
region of 0.1 pc of diameter, b) A^mm vs the non-thermal veloc¬ 
ity dispersion as inferred from VLA NH3(1,1) data, c) Nmm vs 
the non-thermal velocity dispersion as inferred from single-dish 
N 2 H'*" (1-0) data. 


Then, the total (thermal -h non-thermal) velocity disper¬ 
sion is calculated by adding quadratically the thermal and 
non-thermal components, using for the thermal component 
a molecular weight of 2.3, which corresponds to the sound 
speed and thus: criD,tot = Mach number 

M is calculated as Cr 3 D,nth/Cs, with CT3D,nth = v/3criD,nth- 
The resulting Mach numbers range from ~ 2 to 7. 

In Fig. □ we plot the number of millimetre sources 
within a field of view of 0.1 pc in diameter, Nmm (a proxy 
to the fragmentation level), as function of (a) the core den¬ 
sity within a region of 0.1 pc of diameter (modelled in Palau 
et al. 2014) —upper panel—, (b) the non-thermal velocity 
dispersion measured from NH 3 ( 1 , 1 ) (interferometric data) 
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—middle panel—, and from N 2 !!■*■( 1 - 0 ) (single-dish data) 
—lower panel—. The figure shows that, while there is a cor¬ 
relation between A^mm and the density within 0.1 pc (cor¬ 
relation coefficient: 0.89), there is no clear trend between 
A^mm and the velocity dispersion measured with NH3 (cor¬ 
relation coefficient: 0.27), nor with N 2 H''' (correlation coef¬ 
ficient: 0.35). This snggests that the velocity dispersion of 
the massive dense cores might not be a crncial ingredient in 
determining the fragmentation level. 


(Figs. [ 2}3 and [ 2 ^), and including only non-thermal snpport 
but taking into account that large-scale supersonic flows 
compress the gas and generate density enhancements, which 
are the ones that proceed to collapse (Mac Low & Klessen 
2004). In this case, the ‘effective density’ is obtained by mnl- 
tiplying the average density of the core by the square of the 
Mach number: 


1- 71 Tconv.flows 1 
2^ Jeans 

= 0.8255 

ClD.nth 

3 

riHa A4^ 

L Mq 

Lo. 188 kms~M 


. 10 ® cm~®. 


3.2 Fragmentation level vs ‘turbulent’ Jeans 
number 


To further compare the role of the physical properties of 
the cores (density, temperature, velocity dispersion) in de¬ 
termining the fragmentation, we estimated the expected 
number of fragments under different assumptions for the 
gravitationally unstable mass, to which, for convenience, 
we continue referring as a ‘Jeans’ mass in general. To es¬ 
timate the Jeans mass, we started from the Jeans length, 


Ljeans = y (c. g., Kippeuhahn, Weigert & Weiss 2012), 


and assumed spherical symmetry, Mjeans = ^Peff ( ) , 

where Ceff is the ‘effective sound speed’, peff is the ‘effective 
density’, and G is the gravitational constant. Therefore: 


^ _3 .- 1/2 

JWJeans — QQ 3/2 Peff ' 

First, we have searched for a correlation between the 
observed fragmentation level A^mm and the expected num¬ 
ber of fragments in a turbulent support scenario (e. g., Mac 
Low & Klessen 2004). Thus, we have computed the expected 
mass of the fragments by assuming that the critical ‘Jeans’ 
mass is determined by non-thermal (‘turbulent’) support, 
where the ‘effective sound speed’ Ceft corresponds to the 
non-thermal component of the observed velocity dispersion 
ciD.nth. This, in practical units and using the number den¬ 
sity of H 2 molecules (as calculated by the model in Palau 
et al. 2014, and using a molecular weight of 2.8), can be 
written as: 


M? 


L Mq 


0.8255 


fTlD.nth 

3 

31H2 

Lo. 188 kms“% 


. 10® cm“®. 


( 2 ) 


Then, the number of expected fragments, A^jeans, is es¬ 
timated from the ratio of the mass of the core inside a region 
of 0.1 pc of diameter, Mo.iptQ, and the Jeans mass, Mjeans: 


A^Jeans — 


Mo.lpc 

Tfjeans 


(3) 


The result is presented in Figs. [2^ and[2}i (for NH 3 ( 1 , 1 ) 
and N 2 H'^ ( 1 “ 0 ), respectively). 

We also estimated the Jeans mass including both ther¬ 
mal and non-thermal support: 



= 0.8255 

UlD.tot 

3 

31H2 

L Mq 

Lo. 188 kms“% 


. 10 ® cm“®. 


(4) 


The mass inside a region of 0.1 pc of diameter is typically ~ 10% 
of the total mass of the core (given in Table 4 of Palau et al. 2013) 
and is only marginally correlated to the mass of the core, in part 
because far from the central part of the core, the core departs from 
sphericity, a basic assumption of the core modelling of Palau et 
al. (2014). 


(Figs. [2J: and[2f). 

Figs. [2^ (NH 3 ) and[2}l (N 2 H^) reveal a very weak cor¬ 
relation of Nmm with Vjeans (correlation coefficient of 0.24 
and 0.23, respectively), and a slope 1 (0.14 ± 0.18 for 
NH 3 , and 0.22T0.25 for N 2 H''', see Table[2}. For the case of 
non-thermal-|-thermal support the situation is very similar 
(Figs. [213 and (2^), and for the case of ‘density enhanced by 
turbulence’ (equation (5) and Figs. (2}: and (2!), the corre¬ 
lation coefficient increases up to 0.50 and the slope up to 
0.34 ± 0.18 (Table [2J, with respect to the case of turbulence 
providing only support (equations (2) and (4)). This was 
expected, because when taking into account the density en¬ 
hancements produced by turbulence the role of turbulence 
providing support becomes less important and the correla¬ 
tion slightly improves. 

Therefore, we show that considering the non-thermal 
motions as a form of support does not provide a good cor¬ 
relation between the expected number of fragments and the 
observed number in any of the cases considered. More im¬ 
portantly, we note that the turbulent Jeans number for the 
majority of the cores is less than or similar to unity in all 
three cases, which would imply that, if turbulent support 
were active, these cores should not fragment at all, contrary 
to what is observed. As we will see in the next Section, this 
would imply a Core Formation Efficiency (CFE, see below) 
> 100 %, which is meaningless. 


3.3 Fragmentation level vs thermal Jeans number 

Given the poor correlations found between the observed 
number of fragments and the expected number of fragments 
in case of turbulent support, we considered only thermal 
support (no contribution from ‘turbulence’). In this case, 
the ‘effective sound speed’ Ceff simply corresponds to the 
sound speed of the gas, which can be written in terms of the 
kinetic temperature, and the Jeans mass is hnally written 
as: 


M, 


L Mq 


0.6285 


r ^ 

3/2 

riHa 

[iokJ 


. 10 ® cm“®. 


( 6 ) 


This was done using two different assumptions for the 
temperature. First, we used a fixed temperature of ~ 20 K, 
as a first approximation to the ‘initial’ (i.e., before being 
heated by the protostellar feedback) temperature of the 
dense core (e. g., Sanchez-Monge et al. 2013). Second, we 
used the average temperature estimated for each core within 
a region of 0.1 pc of diameter, To.ipc (ranging from 25 to 
150 K, Table[T]). This assumption should give an upper limit 
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Figure 2. ‘Fragmentation level’ {Nmm) vs Jeans number. Top panels: Jeans number calculated using the velocity dispersions estimated 
from NH3, as explained in Section 3.2. Middle panels: idem but using the velocity dispersions estimated from N2H"^ (Section 3.2). 
Bottom panels: Jeans number calculated considering only thermal support, calculated either using a fixed temperature of 20 K for all 
the cores (panels ‘g’ and ‘i’), or the temperature inferred from the core modelling presented in Palau et al. (2014; panel ‘h’; see Section 
3.3). In all panels, the Jeans number is calculated using the average density inside a region of 0.1 pc of diameter (as explained in the 
table notes of Table 1; see also Palau et al. 2014), except for panels on the right, where the average density has been multiplied by the 
square of the Mach number (following Mac Low Sz Klessen 2004). In all panels, the dotted black line represents the one-to-one relation, 
for a core formation efficiency (CFE) of unity. For ‘g’, ‘h’, and ‘i’ panels, the red dashed line corresponds to the fit with slope= 1 used 
to infer the indicated CFE (3—41%). 


to the temperature at the time when fragmentation took 
plac^. 


^ The stellar feedback should affect the density structure on 
larger timescales compared to the timescale when stellar feedback 
modifies the temperature because the first should change through 
mechanical processes while the latter changes through radiative 
processes. In addition, the massive dense cores of our sample are 


in similar evolutionary stages, having not developed UCHII re¬ 
gions yet (see Palau et al. 2014 for a more detailed discussion), 
and the Jeans mass depends more strongly on temperature than 
on density. For these reasons we consider that the density struc¬ 
ture of the massive dense cores in our sample is a reasonable 
approach to the density structure at the time of fragmentation. 
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Table 2. Linear fits to the Afmm vs Afjeans relations of Fig. [2] which correspond to different cases of core support 


Support ^ 

equation ^ 

panels 
Fig. 2 ^ 

correlation 
coefficient ^ 

slope ^ 

CFE == 

(%) 

^2 c 

turbulent (NH 3 ) 

( 2 ) 

a, d 

0.24 

0.14 ±0.18 

> 100 % 

_ 

turbulent (N 2 H+) 

( 2 ) 

a, d 

0.23 

0.22 ±0.25 

> 100 % 

- 

turbulent nAi^ (NH 3 ) 

(5) 

c, f 

0.50 

0.34 ±0.18 

> 100 % 

- 

turbulent nA4^ (N 2 H+) 

(5) 

c, f 

0.36 

0.33 ±0.24 

> 100 % 

- 

thermal T=20K 

( 6 ) 

g 

0.72 

0.60 ±0.14 

13% 

1.27 

thermal varying T 

( 6 ) 

h 

0.57 

0.34 ±0.12 

41% 

3.34 

thermal T=20K nM'^ (NH 3 ) 

( 8 ) 

- 

0.69 

0.35 ±0.11 

3.7% 

2.37 

thermal T=20K nM^ (N 2 H+) 

( 8 ) 

i 

0.64 

0.47 ±0.15 

3.3% 

1.49 


® Fits are performed for the relations of the panels (of Fig. 2) indicated in column (3), which correspond to A^jeans estimated using the 
equations given in column ( 2 ). 

Correlation coefficient and slope of a linear fit with two free parameters. The slope should be close to one if the form of support correctly 
described the observations. 

^ Core Formation Efficiency inferred forcing a linear fit with slope = 1, and the corresponding 


The results are plotted in Fig. [2^,h. Fig. [2^ shows a 
correlation of Afmm and A^jeans, with a slope of 0.60 ± 0.14, 
clearly larger and closer to 1 than the slope obtained for the 
turbulence-supported case (0.2-0.3, Table 2). In this panel, 
the temperature is fixed for all the cores and equal to 20 K. 
The data are clearly offset with respect to the one-to-one 
relation (dotted black line), which can be explained if only 
a percentage of the total mass of the core is converted into 
compact fragments. We dehne the Core Formation Efficiency 
(CFE) as the fraction of mass of a dense core found in (pre- 
and proto-stellar) compact fragments (as in Bontemps et al. 
2010 ), and in this case: 


^ Mo.ipc CEE 

Afjeans 


(7) 


Thus, we fitted a line with slope 1 (dashed red line in Fig. [2]) 
and the offset should be a first approximation to the CFE. 
By doing this for the dataset of Fig. we found a CFE 
of 13%, with a correlation coefficient of 0.72. This value for 
the CEE is fully consistent with the independent direct mea¬ 
surements of the CFE by Bontemps et al. (2010) and Palau 
et al. (2013), who estimated this quantity by dividing the 
total mass in compact fragments (detected with an interfer¬ 
ometer in an extended configuration) by the mass of the core 
(measured with a single-dish). Our inferred CFE is also sim¬ 
ilar to those found by Louvet et al. (2014) in the W43-MM1 
region. 

We additionally estimated A^jeans using the different av¬ 
erage temperatures inferred for each core (Table [T]), and the 
result is shown in Fig.[2]i. In this approach, A^jeans is smaller 
(compared to the previous case of fixed temperature equal to 
20 K), because the Jeans masses are larger due to the higher 
adopted temperatures, and hence the inferred CFE is larger 
as well. The effect of using these higher temperatures is to 
predict too small a number of fragments (too small Aijeans), 
especially for the two extreme cases (cores 13 and 14) which 
are also the most luminous regions. In this case we obtained 
a correlation coefficient of 0.57, and a CFE of 41%. 

Einally, we calculated the Jeans mass considering that 
turbulence is only producing regions of higher density, but 


not providing support against gravity (with the latter being 
only thermal, e.g., Padoan & Nordlund 2002): 


P n^conv.flows —th- 
Jeans 

= 0.6285 

T 

3/2 

■ riHa ' 

- 1/2 

L Mq 

[iokJ 


. 10® cm“® . 



We studied this case using Mach numbers calculated from 
NHa and N2H’'' data (Table 1), and the results are listed 
in Table 2. The fit performed using NH3 (to estimate the 
Mach number) has a correlation coefficient very similar to 
the coefficient obtained for pure thermal support at a fixed 
temperature of 20 K in Fig. 2g, but the slope of the linear 
fit is significantly smaller (0.35 ± 0.11) and thus deviates 
more strongly from the one-to-one relation. As for the fit 
performed using N 2 H'^ (to estimate the Mach number), the 
slope is more similar to the slope in Fig. 2g. In Fig. 2i we 
show the case of N 2 H’^ only for clarity. We also performed a 
linear fit forcing the slope to 1 to infer the CFE in this case, 
which is around 3% for both NH3 and N2H''' (see Table 2 
and Fig. 2i). The CFE is very low because the densities in 
this case are higher (by M^) and the Jeans mass decreases 
resulting in a very high number of expected fragments. 

Overall, the best correlation between A^mm and Ajeans 
is found for the case of pure thermal support adopting a 
temperature of ~ 20 K for all the cores and with no mod¬ 
ification of the density by the Mach number (Eig. 2g). In 
addition, also for this case the slope in the Amm vs Ajeans 
relation is closest to 1 (see Table [2J. 


3.4 Fragment masses 

We have estimated the masses of the fragments identified 
in each massive dense core by assuming the temperature at 
the distance of the fragment (from the core centeiQ) as pro¬ 
vided by our modelled envelopes. We used the dust opacity 
law of Ossenkopf & Henning (1994, icy mantles for densi¬ 
ties ~ 10® cm“®). The results are shown in Eig. [3] About 


® The core center is taken as the peak of the millime¬ 
tre/submillimetre emission observed with a single-dish (see Palau 
et al. 2014 for further details). 
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Figure 3. Histogram of fragments masses. The red solid line 
corresponds to the Jeans mass, averaged over all the sample, 
assuming thermal support at a fixed temperature of 20 K; the 
red dashed line corresponds to the Jeans mass, averaged over all 
the sample, assuming thermal support at the average tempera¬ 
ture measured for each core inside a region of 0.1 pc of diameter 
(which should be an upper limit to the temperature at the time of 
fragmentation); and the green dotted-dashed line corresponds to 
the Jeans mass, averaged over all the sample, assuming turbulent 
support as described in equation (5) (Mac Low &; Klessen 2004). 
For the other prescriptions of turbulent support lSection l3.2l l. we 
obtain Jeans masses up to ~ 80 Mq. 


45 fragments (out of 75) present masses < 1 Mq . Although 
these masses have been inferred from interferometric obser¬ 
vations using extended configurations, and thus part of the 
flux must have been filtered out by the interferometer, we 
estimate that the missed flux is probably not larger than 
a factor of ~ 2 (the reason for this is that any fragment 
detected by the interferometer at such extended configura¬ 
tions must be intrinsically compact). For example, for the 
case of 122198, Sanchez-Monge et al. (2010) report a flux 
density at 1.3 mm using the SMA in compact configuration 
of ~ 500 mjy, while Palau et al. (2013) report a flux density 
at 1.3 mm using the PdBl in its most extended configura¬ 
tion of ~ 270 mJy (adding MM2 and MM2S). Thus, as a 
first reasonable approach, one might conclude that most of 
the fragments in our sample must have masses of the order 
of 1-2 Mq, or below. 

On the other hand, the Jeans masses calculated for each 
region considering only thermal support for a fixed temper¬ 
ature of 20 K (shown as a solid line in Fig. [3| range from 
0.4 to 1 Mq, with an average value of 0.75 Mq. If we use 
our second approach for the temperature, i. e., use the av¬ 
erage temperature estimated for each core inside a region 
of 0.1 pc of diameter (Table [T] without including cores 13 
and 14, whose temperature is clearly affected by stellar feed¬ 
back), we find Jeans masses in the range from 1 to 5 Mq, 
with an average Jeans mass for all the regions of 2.2 Mq 
(shown as a dashed line in Fig. [3|. Since this last case should 
yield an upper limit to the temperature at the time of frag¬ 
mentation, the Jeans mass in this case is also an upper limit. 


For the case of average Jeans masses including the turbu¬ 
lence as a support term (equations 2, 4, and 5), these range 
from 10 Mq (equation 5) up to ~ 80 Mq (equations 2 and 
4). Therefore, we conclude that the typical masses of most 
of the fragments, around 1 Mq, are in good agreement with 
the Jeans mass considering pure thermal support only, with 
no need of additional forms of support. 


4 DISCUSSION AND CONCLUSIONS 

In the present work we have investigated which process, tur¬ 
bulence or gravity, is controlling the fragmentation level of 
massive dense cores at a scale of ~ 0.1 pc. We have used 
a sample of 19 massive cores, previously presented in Palau 
et al. (2013, 2014), to show that the fragmentation of these 
objects seems to be controlled mainly by thermal Jeans frag¬ 
mentation. Specifically, we have shown that the fragmen¬ 
tation level, Amm —measured as the number of compact 
fragments within a core—, presents a significantly better 
correlation with the core density than with its non-thermal 
velocity dispersion. Correspondingly, we have shown that 
Amm correlates linearly with the number of fragments ex¬ 
pected from simple thermal Jeans fragmentation, with a core 
formation efficiency, CFE, around 13%, while, instead, as¬ 
suming turbulence-dominated fragmentation, the majority 
of the cores should not fragment, contrary to what is ob¬ 
served. Finally, we have given hints that the masses of most 
of the fragments in our sample seem to be of the order of the 
Jeans mass calculated considering only pure thermal sup¬ 
port. 


4.1 Comparison to previous works 

It is important to point out that our results are not incon¬ 
sistent with those of Zhang et al. (2009), Pillai et al. (2011), 
Wang et al. (2011, 2014), and Lu et al. (2015). Those au¬ 
thors concluded that the fragments within massive cores of 
infrared-dark clouds have masses significantly larger than 
the thermal Jeans mass, and consistent with the turbulent 
Jeans mass instead. However, in most of these clouds, this 
could be related to a sensitivity and spatial resolution is¬ 
sue due to the large distance of these infrared-dark clouds 
(ranging from 3.3 to 7.4 kpc). For example, in four (out of 
six) of those clouds, the mass sensitivity is > 2 Mq (above 
the Jeans mass) and the spatial resolution is > 5000 au, 
while the massive dense cores studied here are all observed 
with mass sensitivities < 1 Mq, and spatial resolutions 
~ 1000 au. 

The most puzzling clouds are G28.34-1-0.06 PI (Zhang 
et al. 2015), and the Snake (Gll.ll—0.12 P6, Wang et al. 
2014). In these two clouds, observed down to subsolar mass 
sensitivities, there seems to be a lack of low-mass fragments, 
suggesting that the bulk of low-mass fragments have not 
formed yet at such earlier stages, or that the low-mass frag¬ 
ments form outside the core and follow the global collapse 
of the cloud. However, other recent multiwavelength studies 
focused on the stellar content of infrared-dark clouds show 
that most of the protostars formed in these clouds are of low- 
mass (< 2 Mq-, e. g., Samal et al. 2015), and thus this needs 
to be further investigated. Our study, carried out toward a 
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sample of 19 regions more evolved than those in infrared- 
dark clouds, shows that most of the fragments detected in 
our sample are of low mass, and consistent with the ther¬ 
mal Jeans mass, indicating that at these stages the low-mass 
fragments do already exist. If the lack of low-mass fragments 
in infrared-dark clouds is confirmed in future observations, 
our data suggest that the duration of the stage when the low- 
mass fragments are formed is quite short. This is consistent 
with the extremely non-linear nature of the gravitational 
collapse (see, e. g.. Fig. 1 (bottom-left) of Toala, Vazquez- 
Semadeni, & Gomez 2012, and Fig. 1 of Zamora-Aviles & 
Vazquez-Semadeni 2014). 

Finally, while it is possible that, in order to form 
the most massive fragments, additional compression mecha¬ 
nisms besides gravity may be necessary, the bulk of the frag¬ 
mentation process in cores actively forming stars seems to be 
dominated by gravity rather than by turbulence. This is con¬ 
sistent with recent claims that the bulk of the non-thermal 
motions in clouds and cores may be dominated by infall 
rather than by random turbulence (e. g., Vazquez-Semadeni 
et al. 2008, Schneider et al. 2010, Ballesteros-Paredes et al. 
2011a, Peretto et al. 2013, Gonzalez-Samaniego et al. 2014; 
see also the review by Vazquez-Semadeni 2015). 

4.2 Physical implications 

Our results suggest that Amm does not seem to depend sig¬ 
nificantly on the internal supersonic motions of the core, 
and are thus contrary to the widespread notion that support 
against gravity is necessary and that turbulence is the main 
physical process providing it and causing the fragmentation 
of molecular clouds. Since non-thermal supersonic motions 
are indeed observed in massive dense cores, but they do 
not seem to be random enough to act as a pressure against 
gravity, we propose that the observed ‘turbulence’ cannot 
be used to define a ‘turbulent-Jeans’ mass. 

Although turbulence may very well play a crucial role 
in the formation of the seeds of what eventually will grow as 
cores, as demonstrated by the early evolution of molecular 
clouds in numerical simulations (e.g., Clark & Bonnell 2005; 
Vazquez-Semadeni et al. 2007, Heitsch & Hartmann 2008), 
one possible interpretation of our results is that the fragmen¬ 
tation process in star-forming regions is controlled mainly 
by gravitational contraction and the ensuing reduction in 
the thermal Jeans mass as the density increases during the 
collapse. Thus, a possibility is that the non-thermal mo¬ 
tions are dominated by infall, produced by the gravitational 
contraction (e.g., Ballesteros-Paredes et al. 2011a). Indeed, 
analysis of the dense regions in simulations of driven, isother¬ 
mal turbulence, indicate that the overdensities tend to have 
velocity fields with a net negative divergence (i.e. a conver¬ 
gence), rather than being completely random with zero or 
positive net divergence, as would be necessary for the bulk 
motions to exert a ‘turbulent pressure’ capable of opposing 
the self-gravity of the overdensities (Vazquez-Semadeni, et 
al. 2008, Gonzalez-Samaniego et al. 2014). Another possibil¬ 
ity is that the non-thermal motions are strongly affected by 
stellar feedback, but in this case the effect on the clouds may 
be disruptive rather than supportive (Colfn et al. 2013). 

Our work supports the notion that non-thermal mo¬ 
tions cannot be treated as capable of exerting a net turbu¬ 
lent pressure that can provide support against gravity and 


stabilize the cores, since we have found no evidence that 
the ‘turbulent Jeans mass’ plays any significant role in the 
fragmentation of the cores. If this view is correct, then the¬ 
oretical models based on the hypothesis of turbulent sup¬ 
port (e.g., McKee & Tan 2003; Krumholz & McKee 2005) 
should be revised, as well as observational works that over¬ 
simplify the role of turbulence and estimate the Jeans mass 
by using an ‘effective sound speed’ corresponding to the non- 
thermal velocity dispersion. Glearly, a detailed comparison 
with simulations is needed to understand the origin of the 
non-thermal motions in massive dense cores and their role in 
the fragmentation process of molecular clouds. These sim¬ 
ulations would help establish more clearly that gravity is 
indeed controlling fragmentation in massive dense cores at 
0.1 pc scales, and even at scales of the entire molecular cloud 
once the clouds are well developed, as suggested by several 
authors (e.g., Clark & Bonnell 2005; Vazquez-Semadeni et 
al. 2007). 
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APPENDIX A: FULL NAMES AND 
COORDINATES OF THE SAMPLE 

In Table Al we provide the full names or other names, co¬ 
ordinates and distances for the sources in our sample. 


APPENDIX B: NH3 AND N2H+ SPECTRA 

In this appendix, we present the NH3(1,1) and N2H^(l-0) 
spectra used in this work to estimate the non-thermal line 
widths for each massive dense core. The spectra, together 
with the hyperfine fits done using the CLASS program of 
GILDAS, are presented in Figs. [BTI and IB2] and result from 
a compilation of data already published (see Section 2 for 
references) or reduced from the VLA archives. 
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Table Al. Main properties of the sample of massive dense cores studied in this work 


Short 

Full 

Position^ 

Distance 

Tbol*’ 

Mobs'" 

name 

name 

a(J2000) 

(5(J2000) 

(kpc) 

(Lq) 

(Mq) 

1-IC1396N 

IRAS 21391+5802 

21:40:41.71 

+58:16:12.8 

0.75 

290 

78 

2-122198 

IRAS 22198+6336 

22:21:26.78 

+63:51:37.6 

0.76 

340 

115 

3-N2071-1 

NGC2071 

05:47:04.78 

+00:21:43.1 

0.42 

440 

80 

4-N7129-2 

NGC 7129-FIRS2 

21:43:01.68 

+66:03:23.6 

1.25 

460 

81 

5-CB3-mm 

CB3-mm 

00:28:42.70 

+56:42:06.8 

2.50 

700 

169 

6-I22172N 

IRAS 22172+ 5549-N 

22:19:08.60 

+56:05:02.0 

2.40 

830 

119 

7-OMC-lS 

OMC-IS 

05:35:14.00 

-05:24:00.0 

0.45 

2000 

158 

8-A5142 

AFGL5142 

05:30:48.02 

+33:47:54.5 

1.80 

2200 

356 

9-I05358NE 

IRAS 05358+3543-NE 

05:39:13.07 

+35:45:50.5 

1.80 

3100 

1480 

10-120126 

IRAS 20126+4104 

20:14:26.04 

+41:13:32.5 

1.64 

8900 

68 

11-122134 

IRAS 22134+5834 

22:15:09.23 

+58:49:08.9 

2.60 

11800 

222 

12-HH80-81 

IRAS 18162-2048 

18:19:12.10 

-20:47:30.0 

1.70 

21900 

333 

13-W3IRS5 

W3-IRS5 

02:25:40.77 

+62:05:52.5 

1.95 

140000 

971 

14-A2591 

AFGL 2591 

20:29:24.90 

+40:11:19.5 

3.00 

190000 

784 

15-Cyg-N53 

Cygnus X-N53 

20:39:03.10 

+42:25:50.0 

1.40 

300 

675 

16-Cyg-N12 

CygnusX-N12 

20:36:57.40 

+42:11:27.5 

1.40 

320 

622 

17-Cyg-N63 

Cygnus X-N63 

20:40:05.20 

+41:32:12.0 

1.40 

470 

160 

18-Cyg-N48 

Cygnus X-N48 

20:39:01.50 

+42:22:04.0 

1.40 

4400 

865 

19-DR21-OH 

DR21-OH 

20:39:01.00 

+42:22:46.0 

1.40 

10000 

526 


® Approximate position of the center of the field of view (corresponding to a region of ~ 0.1 pc of diameter) where fragmentation was 
assessed in Palau et al. (2013, 2014). 

Bolometric luminosity as given in Table 1 of Palau et al. (2014). 

Mobs is ths mass computed analytically from the model of Palau et al. (2014), integrating until the radius where the density profile 
could be measured for each source. Note that for 122198 we present here an updated version of the model (see Appendix C). 


APPENDIX C: NEW DENSITY AND 
TEMPERATURE DETERMINATION FOR 
IRAS 22198+6336 

In Palau et al. (2014) the original images of 
IRAS 22198+6336 (122198) published by Jenness et al. 
(1995) using the JCMT were not available, and the images 
were digitized. The images are now available and we have 
re-done the fit, with the additional difference (with respect 
to Palau et al. 2014) that the main beams assumed here 
are 7 and 14 arcsec at 450 and 850 /tm, respectively, with 
no consideration of error beams. The results are presented 
in Fig. [Cl] and Table jCTj 
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Figure Bl. VLA NH3(1,1) spectra averaged over a region of ~ 0.1 pc of diameter, where the fragmentation has been assessed. Red 
lines correspond to the CLASS fits to the hyperfine structure, from which the line width given in Table 1 has been inferred. For regions 
122134 and DR21-OH, two velocity components have been used. 


Table Cl. Best-fit parameters to the radial intensity profiles and Spectral Energy Distribution of IRAS 22198-1-6336 (122198), and 
inferred properties (updated after Palau et al. 2014) 




To^ 

PO^ 

rioK^ 

« b 

' max 


V b 

■^0.1 pc 

^0.1 pc^ 

Source 

/ja 

(K) 

(x 10“^^ g cm“®) 

Xr^ (pc) 

(pc) 

(Mq) 

(gcm-2) 

(xlO® cm“®) 

122198 

1.16 ±0.22 

44 ±4 

3.4 ±0.4 

1.75 ±0.03 0.580 0.39 0.22 

0.31 

115 

0.29 

3.6 


^ Free parameter fitted by the model: {3 is the dust emissivity index; Tq and po ^re the temperature and density at the reference radius, 
1000 AU; p is the density power law index; Xr is the reduced x as defined in equation (6) of Palau et al. (2014). 

^ Parameters inferred (not fitted) from the modeling, q is the temperature power-law index, and rioK is the radius of the core where 
the temperature has dropped down to ~ 10 K; rmax is the radius at the assumed ‘ambient’ density of 5000 cm~^; Mobs is the mass 
computed analytically from the model integrating until the radius where the density profile could be measured for each source. So.ipc 
and no.ipc are the surface density and density inside a region of 0.1 pc of diameter. 
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Figure B2. Single-dish N2H'^(l-0) spectra averaged over a region of 0.1 pc of diameter, where the fragmentation has been assessed. 
Blue lines correspond to the CLASS fits to the hyperfine structure, from which the line width given in Table 1 has been inferred. 







































14 Palau et al. 


r (au) 



Frequency (GHz) 

10 ^ 10 '* 1000 100 10 1 



r (au) 



Frequency (GHz) 
1000 100 


22198 SED fit residuals 


100 1000 10 ^ 
Wavelength (//m) 


Figure Cl. New fit (after Palau et al. 2014) to the radial intensity profiles and Spectral Energy Distribution of IRAS 22198+6336 
(122198) after using original JCMT data of Jenness et al. (1995) work. 












































